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O ' 1- Introduction 

00 ' 

Thermal radiation from surfaces of several radio pulsars has been detected 
recently in soft X-rays by ROSAT and ASCA (see Becker and Trumper, 1998). A 
number of the point-like X-ray sources has also been discovered and identified as 
radio-silent isolated cooling neutron stars (NSs) (see Caraveo et al., 1996; Walter 
et al., 1996; Vasisht et al., 1997). In several cases, the observations seem to be 
better explained if the NSs possess an envelope of matter of low atomic weight, 
presumably hydrogen (Page, 1997), and if NS atmosphere models are applied to 
fit the data (Pavlov and Zavlin, 1998). There are indications that this may be the 
case for NSs with different magnetic fields B: from weak, B -C 10 10 G (Rajagopal 
and Romani, 1996; Zavlin and Pavlov, 1998; Zavlin et al., 1998), to strong, B ~ 
10 n ~ 13 G (Page et al., 1996; Zavlin et al., 1998), and superstrong, B > 10 14 G 
(Heyl and Hernquist, 1997). To justfy this, one needs in further observations 
and advanced atmosphere models of NSs with high B. Thermal motion of atoms 
accross the magnetic field induces an electric field in the frame of the atom. This 
affects atomic structure, atmospheric thermodynamics and opacities (Ventura et 
al., 1992; Pavlov and Meszaros, 1993), being a critical point in the construction of 
models and data interpretation. Here we study this problem for hydrogen plasma 
at temperatures T ~ 10 5 ' 5 — 10 6 K, densities p ~ 10~ 3 — 10 1 g cm -3 , and magnetic 
fields B ~ 10 12 — 10 13 G, typical for atmospheres of middle-aged cooling NSs. We 
construct an analytic model of the plasma free energy and derive a generalized 
Saha equation which is used to obtain the opacities. 
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2. Method and results 

We use a generalization of the method applied previously to the non- 
magnetic case (Potekhin, 1996a). The Helmholtz free energy of a volume V of 
the plasma is F(V, T, {N a }) = + F ; ( d e) + F i ( d ff) + F cx , where the subscripts "id" 
and "ex" refer to the ideal and excess part of F, (e, p, H) relate to the electrons, 
protons, and atoms, respectively, and {N a } is the set of numbers of the particles. 

For protons, (3F^/N P = \n(2n a 2 m X p n p ) + In [l - e _/9Rw< »] - 1, where (3 = 
l/k B T, n a = N a /V, X a = (2TT/3h 2 /m a ) 1 / 2 , a m = (hc/eB) 1 / 2 , u cp = eB/m p c. 
Being the same for free protons and atoms, proton spin and zero-point energy do 
not affect thermodynamics, therefore the corresponding terms are suppressed. 

For electrons, we take into account partial degeneracy and relativity: 
Ffo /V = fm e — P[d\ where the pressure P$ is easily expressed through fitting 
formulae by Blinnikov et al. (1996), and n e is related to the chemical potential 
through a fitting formula by Potekhin (1996b). For F^ p \ we use an analytic ex- 
pression (Chabrier and Potekhin, 1998) based on calculations by Chabrier (1990) 
at B — 0. The latter approximation is justified for the classical plasma (accord- 
ing to the Bohr- van Leeuwen theorem) but fails in the region where B suppresses 
strong electron degeneracy (that occurs only in subphotospheric layers). 

For atoms, (3F^ } = £ s „ / d 2 K ± N su (K ± ) {In [n H X 3 H w sl/ (K ± )/Z w ] - 1}, s 
and v relate to electronic excitations, N su = (\ H /2nh) 2 N H w su e /3x3V /Z w are the 
atomic occupancies per unit phase space of the transverse component K± of the 
pseudomomentum K which characterizes the atomic motion in the magnetic field, 
w S v(Kj_) and Xsv{K±) are the occupation probabilities and binding energies of the 
moving atom, and Z w = (\h/^K) 2 J2su I d 2 Kj_w su (K]_) exp[flxsi>{K±)} is the in- 
ternal partition function. Large values of correspond to the decentered states 
with large dipole moment created by the induced electric field. These states have 
been neglected in previous studies (e.g., Lai and Salpeter, 1997) but our results 
show that they are important at T > 10 5,5 K. Generalization of the unbinding 
free energy (Potekhin, 1996a) to the magnetic case gives (3F^ = (N H + N e ) 
Esul d 2 K ± N su v sl/ /V, where v su (K ± ) = (4tt /3)(rfj>\K ± )) 3 is an effective inter- 
action volume and r ep the root-mean-square electron-proton distance in the atom. 

Minimization of F leads to the generalized Saha equation 

(2na 2 ) 2 

n H = \p\ e — -j 2 — n e n p [1 - exp(-phu cp )] Z^exp(A deg ), (1) 

where \nw su w P(dF&)/dN p + dF&>/dN e - (n e + n H )v su (K ± ), and A dcg = /?[// + 
dfj,/dlnn e — dP^ / dn e ]— ln(2im e X e a 2 n ) is the correction for electron degeneracy. In 
numerical solutions of Eq. (1) we use fitting formulae for \sv and r'ff 1 (Potekhin, 
1998). The equation of state (EOS) is calculated from P = -{dF/dV) T ,{N a }- 
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Fig. 1. Fraction of atoms (solid lines — total, dashed lines — ground-state) that 
contrubute to the opacities at different B shown at the panels. Numbers near curves 
are lgT. At high p neutrals dissociate due to the pressure ionization effect. 
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Fig. 2. Opacities of hydrogen plasma at B = 2.35 x 10 12 G, T = 3.16 x 10 5 K, 

and p = 0.1 g cm" 3 for the radiation with left (left panel), longitudinal (middle 

panel) and right (right panel) polarization. Dashed lines represent the contribution 

of ground-state atoms, dot-dashed lines — the opacity of the charged component, 

and solid line — the sum of both. 

Figure 1 shows the neutral fraction / (uh + n p ) that contributes to the 
atomic opacities (n°ff(K±) ~ n sv (K±) exp[— 4 3 n p v su (K±)}, cf. Potekhin, 1996a). 
The total n H contributing to the EOS can be much larger in the low and high 
density limits. The different dependence of n^ 1 on B at high and low T is ex- 
plained by the competition of change in Xsv and the phase space occupied by 
electrons and atoms with B. The opacities shown in Fig. 2 are obtained using the 
bound-bound and bound-free cross sections calculated previously (Potekhin and 
Pavlov, 1997 and ref. therein), and absorption coefficient of the ionized fraction 
calculated following Shibanov et al. (1992). The moving bound species contribute 
largely in the opacities of typical cooling NS atmospheres, their spectra are con- 
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siderably different from those of the static-atom approximation (cf. Pavlov et al., 
1995). This impels us to revise the previous atmosphere models at T < 10 6 K. To 
do that properly one has to include the absorption by excited atoms, which can 
be substantial at relatively high T ~ 10 6 K (see Fig. 1), and the bound- quasifree 
transitions (cf. Stehle and Jacuemot, 1993) which are expected to fill the gap seen 
in the absorption for circular polarizations between b-b and b-f contributions. 
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DFG-RFBR grant 96-02-00177G, and INTAS grants 94-3834 and 96-0542. YS 
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